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Combs are a simple caricature of various types of natural branched structures, which belong to 
the category of loopless graphs and consist of a backbone and branches. We study continuous time 
random walks on combs and present a generic method to obtain their transport properties. The 
random walk along the branches may be biased, and we account for the effect of the branches by 
renormalizing the waiting time probability distribution function for the motion along the backbone. 
We analyze the overall diffusion properties along the backbone and hnd normal diffusion, anomalous 
diffusion, and stochastic localization (diffusion failure), respectively, depending on the characteristics 
of the continuous time random walk along the branches. 

PACS numbers: 05.40.-a 


I. INTRODUCTION 

Random walks often provide the underlying meso¬ 
scopic mechanism for tran^ort phenomena in physics, 
chemistry and biology [ll-l3|. A wide class of random 
walks give rise to normal diffusion, where the mean- 
square displacement (MSD), ((Ar)^(t)), grows linearly 
with time t for long times. In many important appli¬ 
cations, however, the MSD behaves like ((A?/^(t)) oc M, 
with 7 1, and the diffusion is anomalous [llQ- Anoma¬ 

lous diffusion can be modelled by various classes of ran¬ 
dom walks . We focus on the important class of contin¬ 
uous time random walks (CTRWs) [U, A specific fea¬ 
ture of a CTRW is that a walker waits for a random time 
T between any two successive jumps. These waiting times 
are random independent variables with a probability dis¬ 
tribution function (PDF) /(r), and the tail of the PDF 
determines if the transport is diffusive (7 = 1) or sub- 
diffusive (7 < 1). Heavy-tailed waiting time PDFs give 
rise to subdiffusion. Realistic models of the waiting time 
PDF have been formulated for transport in disordered 
materials with fractal and ramified architecture, such as 
porous discrete media @ and comb and dendritic poly¬ 
mers and for transport in crowded environments 

i- 

A simple caricature of various types of natural 
branched structures that belong to the category of loop¬ 
less graphs is a comb model (see Fig.[T]). The comb model 
was introduced to understand anomalous transport in 
percolation clusters [lol - [T^ . Now, comb-like models are 
widely employed to describe various experimental appli¬ 
cations. These models have proven useful to describe 
the transport along spiny dendrites [3 [3, percolation 
clusters with dangling bonds 3, diffusion of drugs in 
the circulatory system [T^ , energy transfer in comb poly¬ 
mers [1,0] and dendritic polymers Q , diffusion in porous 
materials [3-[3|) the influence of vegetation architec¬ 
ture on the diffusion of insects on plant surfaces 
and many other interdisciplinary applications. Random 


walks on comb structures provide a geometrical explana¬ 
tion of anomalous diffusion. 

More general combs have been studied recently. For 
example, a numerical study of the encounter problem 
of two walkers in branched structures shows that the 
topological het erog eneity of the structure can play an 
important role [20j. Another example is the occupation 
time statistics for random walkers on combs where the 
branches can be regarded as independent complex struc¬ 
tures, namely fractal or other ramified branches [2l|. Fi¬ 
nally, we want to mention studies to understand the dif¬ 
fusion mechanism along a variety of branched systems, 
where scaling arguments, verified by numerical simula¬ 
tions, have been able to predict how the MSD grows with 
time p^ . 

Diffusion on comb structures has also been studied by 
macroscopic approaches, based on Fokker-Planck equa¬ 
tions 3: which have been applied to describe diffu¬ 
sive properties in discrete systems, such as porous dis¬ 
crete media Q, infiltration of diffusing particles from 
one material into another [2^ , and superdiffusion due to 
the presence of inhomogeneous convection flow [3 . [3| . 
Other macroscopic descriptions, based on renormalizing 
the waiting time PDF for jumps along the backbone 
to take into account the transport along the branches 
[3| . have been found useful to model continuous-time- 
reaction-transport processes (3| and human migrations 
along river networks [28l| . 

Kahng and Redner provided a mesoscopic, probabilis¬ 
tic description of random walks on combs, by using the 
successive decimation of the discrete-time Master equa¬ 
tion to obtain a mesoscopic balance equations for the 
probability of the walker to be at a given node at a given 
time [2^. A mesoscopic description is necessary for an ac¬ 
curate description of the transport properties, such as the 
diffusion coefficient or the mean visiting time in a branch, 
in terms of the microscopic parameters that characterize 
the random walk process. 

Here we obtain transport quantities within the frame- 
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work of the CTRW formalism. We assume that the 
motion along the backbone and the branches is non- 
Markovian and that the motion along the branches can 
be non-isotropic. We reduce the dynamic effect of the 
branches to a waiting time PDF for the motion along the 
backbone by using the decimation method of Kahng and 
Redner. The time spent by the walker between its en¬ 
try into the branches and its return to the backbone for 
the first time is treated as a contribution to the effective 
waiting time at the node where the branch crosses the 
backbone. 

The paper is organized as follows. In Sec.|TT]we formu¬ 
late the mesoscopic description of the random walk on 
the comb and reduce walker’s motion to an effective mo¬ 
tion along the backbone only with a renormalized waiting 
time PDF for the backbone nodes. Sec. IIIII deals with the 
MSD of the effective backbone motion, derives the effec¬ 
tive diffusion coefficient, and establishes the conditions 
for normal diffusion, anomalous diffusion, and stochastic 
localization (diffusion failure) in terms of the num¬ 
ber of branch nodes and the degree of bias of the motion 
along the branches. We provide details of the numerical 
calculations in Sec. ?? and summarize our results and 
discuss their implications in Sec. El 


II. MESOSCOPIC DESCRIPTION 


The simplest comb model, shown in Fig. [T1 is formed 
by a principal axis, called the backbone, which is a one¬ 
dimensional lattice with spacing a, and identical branches 
that cross the backbone perpendicularly at each node. 
The walker moves through the comb by performing jumps 








backbone 







branch 

a 


FIG. 1: Comb structure consisting of a backbone and 
branches. Each point represents a node where the walker 
may jnmp or wait a random time 

between nearest-neighbor nodes along the backbone or 
along the branches. We assume that the walker performs 
isotropic jumps along the backbone, but the jumps along 
the branches may be biased, for example by an external 
field M- 

We derive the balance equation for the PDF P{x,t) 
of finding the walker at node x on the backbone at time 


t. When the walker arrives at a node, it waits a ran¬ 
dom time T before performing a new jump to the nearest 
node. We assume that the comb is homogeneous, and the 
waiting time PDF at any given node is given by (''')• 
When the walker enters a branch, it spends some time 
moving inside the branch before returning to the back¬ 
bone. This sojourn time can be used to determine an 
effective waiting time PDF (j){T) for the walker’s motion 
along the backbone. In other words, the motion of the 
walker on the comb can be reduced to the effective mo¬ 
tion along a one-dimensional lattice, corresponding to the 
backbone only. This motion is non-Markovian and can 
be described mesoscopically by the Generalized Master 
Equation (GME) 


dP 


f K{t — t')dt'x 
Jo 



x' ,t')^{x')dx' 


P{x,t') 


( 2 . 1 ) 


where K(t) is the memory kernel related to the waiting 
time PDF via its Laplace transform, K{s) = s(/)(s)/[l — 
(^(s)], where s is the Laplace variable. The dispersal ker¬ 
nel $(a;) represents the probability for the walker of per¬ 
forming a jump of length x. If the walker moves isotrop¬ 
ically between nearest neighbors in a one-dimensional 
lattice of spacing a, the dispersal kernel reads $(a;) = 
S{x — a)/2 + 6{x + a)/2. We assume that the walker is 
initially located at a; = 0, i.e., P{x, 0) = (5a;^o with x = ia 
and i = 0, ±1, ±2,..., where Sx,o is the Kronecker delta. 
Then the Laplace transform of the GME ior x ^ 0 reads 

P{x, s) = [P{x - a, s) -I- P{x + a, s)] . (2.2) 

The mesoscopic balance equation for the walker on the 
comb being at node x = ia of the backbone is 


P{x, s) = [p(a; — a,s) + P{x + a, s)] 

+ (1 - 9)0o(s) [P{y = a, s) + P{y = -a, s)]. (2.3) 

Here P{x,s), P{x — a, s), and P{x + a,s) is short¬ 
hand for P(x,y = 0, s), P(x — a, y = 0, s), and 
P{x + a,y = 0,s), and P{y = a,s) and P{y = —a,s) 
stands for P{x,y = a,s) and P{x,y = —a,s). The 
term 0o(s) [P{x — a,s) + P{x + a, s)] /4 corresponds to 
the contribution of the walker arriving at node x = ia 
from the left or from the right with probability 1/4 after 
waiting a random time r with PDF (jjQ (r) at nodes x + a 
or X — a. As shown in Fig. [51 the walker located at the 
jth node of the backbone may jump to the right, left, 
up or down with probability 1/4. We assume that the 
walker moves forward (away from the backbone) along 
the branches with probability q and back to the back¬ 
bone with probability 1 — q. The term 

(1 - g)(/o(s) [Piy = a,s) + P{y = -a, s)] (2.4) 
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in (lO) corresponds the contribution of the walker arriv¬ 
ing at the backbone node x from the first node of the 
upper or lower branch after waiting there a random time 
r with PDF 4>o{t). 

Consider the motion along the upper branches. The 
lower branch dynamics is the same due to the symmetry 
of the comb. The mesoscopic balance equation for the 
first node of the upper branches reads 

P{y = a,s) = s) + (Ao(s)(l - q)P{y = 2a, s). 

(2.5) 

The first term (/)o(s)P(a:, s)/4 corresponds to the contri¬ 
bution of the walker arriving from the backbone, while 
d’o(s)(l — q)P{y = 2a, s) is the contribution of the walker 
jumping from the upper node y = 2ato y = a with prob¬ 
ability 1 — q after waiting a random time r with PDF 
Analogously, we have for the lower branches 

P{V = -a, s) = s)+(j)o{s){l-q)P{y = -2a, s). 

( 2 . 6 ) 

Generalizing (12.51) to any node of the branches located 
between 2a< y < {N — 2)a, we obtain the balance equa¬ 
tion for the upper branches 


P{V, s) = (l)o{s) [qP{y -a,s) + {l- q)P{y + a, s)]. 

(2.7) 

To determine the Laplace transform (j)(s) of the ef¬ 
fective backbone node waiting time PDF, we need to 
determine P{y = a, s) and P{y = —a, s) in (12.3|) in 
terms of P{x,t), so that (12.311 can be cast in the form 
of (j2.2|) . Given (j2.5|) and (12.6|) , this goal can be achieved 
if P{y = 2a, s) and P{y = —2a, s) can be related to 
P{y = a,s) and P{y = —a,s). We proceed as follows. 
The solution of (1^ reads 


P(y,s) = 


where 


A± 


1 ± V^l - 4g(l - g)0§(^ 
2(1 - q)(j)o{s) 


( 2 . 8 ) 


(2.9) 


To find expressions for the quantities Ai and A 2 , whose 
dependence on x and s is not displayed, we apply (12.811 
to the node y = 2a: 


P{y = 2a,s) = AiXl+A 2 Xi. (2.10) 
On the other hand, setting y = 2a in (12.7p . we find 


P{y = 2a, s) = (l)o{s) [qP{y = a, s) 

+ {1 - q)(j)o{s)P{y = 3a,s)], (2.11) 


or 


P{y = 2a, s) - (j)o{s)qP{y = a,s) = 

(j)o{s){l - q)(l)ois)P{y = 3a,s). (2.12) 


Setting y = 3a in ()2.8I) we obtain 


P{y = 2a, s) - q(l)o{s)P{y = a, s) 

= Ms)il-q)[AiXl+A2X^]. (2.13) 

Solving the system of equations (12.1011 and (12.1311 for the 
quantities Ai and A 2 , we obtain 


Ai 


Pjy = 2a, s) - q(j)o{s)P{y = a, s) 
A^ (A+ - A_)((io(s)(l - q) 


X-P{y = 2a, s) 
A^(A+-A_) ’ 


(2.14) 


A2 


-Pjy = 2a, s) + q(j)o{s)P{y = a, s) 

X'i (A+ - A_)(/)o(s)(l - q) 

A-^P(y = 2a, s) 

^ A2_(A+-A_) • 


(2.15) 


A special situation occurs at the end of the branches, 
where we have to impose reflecting boundary conditions, 
i.e.. 


P{y = Na, s) = q(j)o{s)P{y = {N - l)a, s). (2.16) 

The node a.t y = {N — l)a also needs a special balance 
equation (see Fig. [^l. 



FIG. 2: Schematic representation of the possible jumps of a 
walker with the corresponding probabilities. 
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P{y = {N - l)a, s) = q(j)o{s)P{y = {N - 2)a, s) 

+ Ms)P{y = Na,s). (2.17) 

Substituting y = [N — 2)a into (12.71) and considering 
p.l6[l . we can write 

P{y = {N-2)a, s) = h{(j)o{s))P{y = (N-3)a, s), (2.18) 
where 


HMs)) 


#o(s) [1 - #o(s)] 
1 + q(q - 2)(j)l{s) 


(2.19) 


Substituting the solutions from (12.811 . (I2.14L and (12.1511 
into (12.181) . we find 


P{y = 2a, s) = G{q, (j)o{s))P{y = a, s), (2.20) 


where 


G{q,(j)o{s)) 


H{q,(j)o{s)) 


2 #o(s) 


1 + 


1 + H{q,(j)o{s)) 
1 - i7((?,(/)o(s)) 


V'1-49(1 


\X+) X+- h{(j3ois))' 


9)</'o(s) 

( 2 . 21 ) 

( 2 . 22 ) 


For the lower branch we obtain in a similar manner, 


P{y = -2a, s) = G{q, (l)o{s))P{y = -a, s). (2.23) 


We have achieved our goal of expressing P{y = 2a, s) 
and P{y = —2a, s) in terms of P{y = a, s) and P{y = 
—a, s). Substituting (12.201) and (12.2311 into (12.51) and (12.61) 
and using the resulting expressions in (12.311 . we obtain an 
equation of the form (12.21) with 




_ Ms) _ 

(4-g)<('o(g) 

1 - (1 - q)(j)o{s)G{q,(j)o{s)) 


(2.24) 


The Laplace inversion of (12.241) yields (j>{T), which in¬ 
corporates the dynamics along the branches and can be 
understood as the effective waiting time PDF for a walker 
moving along the backbone only. 


III. STATISTICAL PROPERTIES 
A. N finite 

If the local waiting time PDF (j)o{T) has finite moments, 
its Laplace transform reads 0, 0o(s) — 1 — st, in the 
large time limit s —>■ 0, where t is the local mean waiting 
time at each node. Taking the limit s —>■ 0 in (I2.24L we 
obtain the waiting time PDF for the effective backbone 
dynamics, 




FIG. 3: Dimensionless mean waiting time of the effective 
backbone dynamics. 


The mean waiting time (t) is given by 

(t) = ^^[2(l-g)i-^g^ + 4g-3]. (3.2) 

In Fig. [21 we plot the effective mean waiting at a node 
of the backbone dynamics versus N and q. It shows that 
the mean waiting time (t) is a monotonically increasing 
function of both q and N. If the random walk inside 
the branches is isotropic, q = 1/2, one obtains by the 
L’Hopital’s rule from (j3.2ll 

lim (t) = (1 -f 2N) t. (3-3) 

q^l/2 


To determine the diffusion coefficient D for diffusion 
through the comb, we first calculate the MSD. Perform¬ 
ing the Fourier-Laplace transform on m, we obtain 


P{k,s) = 


1 - M) 


s[l — <i)(A:)(/(s)] 

The MSD in Laplace space reads (see, e.g., i) 

<pP{k,s) 


= — lim 
' ' k^O 




(3.4) 


(3.5) 


As mentioned in Sec. m we assume that the motion on 
the backbone is unbiased and that the walker only jumps 
to nearest neighbors. This implies that the kernel $(a;) 
is given <l>(a:) = S{x — a)/2 + S{x + a)/2, and we obtain 
from (13.51) . 




(3.1) 




s [(/(s)-i - 1] 


(3.6) 
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in the large time limit. If the waiting time PDF (j>(t) 
possesses a finite first moment, implies that the 

MSD along the backbone corresponds to normal diffu¬ 
sion = 2Dt. The diffusion coefficient is given by 


D = 



a? 2q — 1 

-f 4q-3' 


(3.7) 


In Fig. m we compare the results provided by (13.7p with 
numerical simulations. 



FIG. 4: Plot of the diffusion coefficient through the comb for 
N = 2, N = ?), and N = Q versus q. Solid curves corre¬ 
spond to exact analytical results given by (pTl) . Numerical 
simulations results are depicted with symbols. 


As Fig. [3] demonstrates, (f) increases monotonically 
with N ior q < 1/2 and saturates at (4g —3)/(2g— 1) for 
N ^ oo. Consequently, the mean waiting time {t) is fi¬ 
nite for N ^ oo; the overall diffusion along the backbone 
is normal. However, for q > 1/2, the mean waiting time 
{t) increases without bound as N increases, and anoma¬ 
lous transport is expected for —>• oo. In Fig. [51 we plot 



t 

FIG. 5: MSD over 2D for N fixed and three different values 
olq- 0.1, 0.25, 05. 0.75 


the result given by (Irfl) for the MSD. The transport is 
diffusive for finite N, regardless q and the specific form 
of (pqiyT), as long as it has finite moments. 

We consider now a an effective waiting time PDF with 
the large-time limit 4‘q{T) ^ with Laplace trans¬ 

form (/o(s) ~ 1 — (stq)'’' and 0 < 7 < 1 , which does not 
possess finite moments. Here tq is a parameter with units 
of time. In this case, the waiting time PDF for the back¬ 
bone dynamics is obtained by simply replacing st with 
(sTo)'’', i.e, (j){s) ~ [l-l-(sTo)"’' (t) /to]“^. Substituting this 
result into (13.61) . we hnd 




a'^To 

it) r(i + 7)’ 


(3.8) 


for large t, where (t) is given by ()3.2I1 . with tq instead of t. 
If the waiting time PDF />o(t) at each node of the comb 
has a power-law tail, then the overall transport along the 
backbone is anomalous. 


B. A —>■ 00 


If the number of nodes of the branches goes to inhnity, 
the mean time spent by the walker visiting a branch in¬ 
creases monotonically, see (USD. However, this does not 
always results in anomalous transport along the overall 
structure as we show below. 

For A —>• 00 , the quotient (A_/A+)^ —>■ 0 and also 
H —>■ 0. We obtain from (12.2111 . 


G{q,(t)o{s)) 


_2#o(5)_ 

1 -b y/1 - 4g(l - 9 )</o(s) 


2#o(5) 

1 + 5(9)’ 

(3.9) 


where we define g{q) = \/l — 4^(1 — q)(/)Q(s) for conve¬ 
nience. Equation (12.241) reduces to 


,, X _ Ms) [1 + 5 ( 9 ) - 2 g(I - q)Ms)] 

2 - (1 + 3<7 - M)Ms) + [2 - (1 - q)Ms)]9iq)' 

(3.10) 

We take the limit s —>■ 0 and consider first the case where 
(j)o{T) has finite moments. Then /)o(s) ~ 1 — st, as s —>■ 0. 
The square root g{q) in (I3.10|l reads 


9 ( 9 ) 


l-2q-^^^si, q<l/2, 

< v^-^(st)3/^, 9 = 1/2, 

-1 + 2,+ tt^st-, ,>1/2, 
2,-1 


(3.11) 


the MSD scaled by the diffusion coefficient. It illustrates 
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and (I3.10p implies that the waiting time PDF is given by 

-1 


(^(s) ~ < 


, 4g-3 - 

2q - 1 


3 ( 7-1 Aq^-Sq+l 


q < 1/2, 

9=1/2, 

-1 

si] , q>l/2. 


q {2q-l)q 

Substituting (13.121) into (13.61) . we find for large t 

2 2 g — 11 


(3.12) 


(x'^it)) = < 


4g — 3 t 
iri 

2 q 


q < 1/2, 

(7 = 1/2, (3.13) 


2 ( 7-1 

where the rate of saturation is 


(l-e-“*), (7 > 1 / 2 , 


0 < 7 < 1, ass—>-0. The MSD in this case can be 
obtained straightforwardly by replacing si with (stq)’^ in 
(13.121) . For large times it reads 


2q — 1 f t 


{x\t)) = 


F(1 + 7 ) dg - 3 Vro^ 

.2 

T-oJ 

Kt/To), 


V2F(1+7/2) \toJ 
a^q{2q- 1 ) 

Tq (V - 3g + 1) 


where 


q < 1/2, 


(7 = 1 / 2 , 


<7 > 1/2, 

(3.15) 


Kt/To) = {t/royE^^j+i 


ftV {2q-ir - 

\toJ V_3(7+1_ 


(3.16) 

is expressed in terms of the the generalized Mittag-Leffler 
function Ea,p{z). We use the following property of inte¬ 
gration of the Mittag-Leffler function [3]] |. 


(29-1)^ 

( 4 ( 7 ^ -3q + l)t' 


(3.14) 


In Fig| 6 ]we compare these results with numerical simu¬ 
lations for N = 10^. For q = 1/2, we obtain the well 
known result of subdiffusive transport with the MSD 
^ y/i. However, for q ^ 1/2, the side branches expe¬ 
rience advection, and the transport is remarkably differ¬ 
ent. Namely, for g > 1/2 the advection is away from 
the backbone along the branches, y —>■ ± 00 . The walker 
is effectively trapped inside the branches, and stochas¬ 
tic localization (diffusion failure) occurs, (a;^(oo)) < cx), 
[^ . For q < 1/2, the advection is towards the backbone. 
It enhances the backbone dynamics and normal diffusion 
takes place. 



FIG. 6 : MSD for three values of q, displaying three different 
behaviors. Solid curves correspond to the results given by 
(13.131) . Symbols are the results of numerical simulations with 
N = lOT 

Consider now the case where the local waiting time 
PDF is (/o(t) ~ i.e., (/o(s) — 1 — (stq)’^ with 


( 6 z“) z^-^dz = (be). (3.17) 

Subdiffusion in the branches results in backbone subdif¬ 
fusion for q < 1/2. For advection away from the back¬ 
bone, q > 1/2, we again find stochastic localization. For 
t/ro >> 1, Ea^p{—at°‘) ^ t““/r(/3 — a) [S^, and conse¬ 
quently y,{t/To) approaches a finite value as t —>■ 00 . 



IV. CONCLUSION 

We have developed a mesoscopic equation for a ran¬ 
dom walk on a regular comb structure given by (lO) and 
(12.241) . The random walk along the branches consists of, 
possibly biased, jumps to the nearest node, while waiting 
at each node for a random time r distributed according 
to the PDF (j)o{T) before proceeding with the next jump. 
The overall dynamics along the branches has been re¬ 
duced to an effective waiting time PDF, given by (I2.24L 
for motion solely along the backbone. We have obtained 
statistical properties, such as the effective mean waiting 
time, (t) for the backbone nodes, and the diffusion coef¬ 
ficient, D, of the overall structure for the case where the 
number of nodes N of the branches is finite or infinite. If 
N is finite and 4>o(x) has finite moments, both (t) and D 
are derived exactly in terms of the bias probability g, the 
number of nodes N on the branch, and the mean waiting 
time probability at each node. In this case the transport 
is always normal diffusion. If 4>q{T) ^ for large 

time, it does not posses finite moments and the MSD of 
the random walker behaves like . If IV is infinite, the 
value of g is decisive. If />o(''‘) has finite moments, the 
diffusion regime is normal if g < 1/2, while the MSD 
behaves like for g = 1/2. If g > 1/2, the MSD 
approaches a constant finite value for large time, corre¬ 
sponding to stochastic localization (diffusion failure). If 
~ T~^~^ for large time, the MSD behaves like 
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for q < 1/2 and like for q = 1/2 Again, stochastic 
localization occurs for q > 1/2. In summary, if the bias 
probability of moving away from the backbone is g > 1/2, 
then stochastic localization occurs, regardless of the other 
characteristic parameters related to the random walk on 
the branches. 
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